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Abstract 


An  overview  of  high-explosive  (HE)  blast  phenomenology,  its  models,  methods,  and  the 
scaling  of  results  are  presented.  Experimental  data  of  HE  blast  damage  to  various  types  of 
structures  are  also  presented,  together  with  some  of  the  HE  blast  structural  damage  models. 
Empirical  vulnerability  models  of  weapon  systems  to  HE  blast  are  also  discussed. 
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1.  Introduction 


To  explain  what  is  meant  by  the  term  “high-explosive  (HE)  blast  damage  mechanism,”  we 
first  need  to  define  the  meaning  we  attach  to  the  words  “blast  wave.”  A  blast  wave  is  a  pressure 
wave  generated  in  the  air  by  the  rapid  release  of  energy  stored  in  some  source  into  the 
surrounding  medium.  There  are  many  sources  of  stored  energy  that  give  rise  to  such  pressure 
waves.  The  stored  energy  in  a  compressed  gas  or  vapor,  either  hot  or  cold,  can  be  such  a  source. 
The  failure  of  a  high-pressure  gas  storage  vessel  or  boiler,  or  the  muzzle  blast  from  a  gun,  is  an 
example  of  explosions  that  give  rise  to  a  blast  wave  in  air.  However,  the  more  usual  sources  of 
energy  for  explosions  and  generation  of  blast  waves  are  either  chemical  or  nuclear  materials  that 
are  capable  of  violent  reactions  when  properly  initiated.  Here  we  shall  confine  our  interest  to 
blast  waves  generated  by  the  detonation  of  HE  charges  (i.e.,  chemical  energy  sources). 

In  order  to  facilitate  the  study  of  HE  blast  waves,  an  explosive  material  is  chosen  as  a 
reference  material  for  all  explosives,  and  the  properties  of  blast  waves  from  spherical  or 
cylindrical  charges  of  this  reference  explosive  are  documented  and  considered  a  standard 
explosive.  Two  explosives  have  been  used  as  such  a  standard,  TNT  and  Pentolite,  and  their 
properties  have  been  studied  and  documented.  The  blast  properties  of  all  other  explosives  are 
usually  derived  by  some  equivalency  to  one  of  the  standard  explosives.  The  blast  properties  of 
TNT  were  documented  by  Swisdak  [1],  mostly  for  use  in  U.S.  Navy  (USN)  publications.  The 
detonation  and  blast  properties  of  Pentolite  were  studied  by  many  researchers  and  documented  by 
Goodman  [2]  for  the  U.S.  Army  Ballistic  Research  Laboratory  (BRL).  Pentolite  is  widely  used 
as  a  standard  explosive  because  of  its  more  reproducible  properties.  The  detonation  and  blast 
properties  of  24  selected  explosives  were  calculated  by  Shear  and  Arbuckle  [3].  Many  other 
publications  exist  on  the  subject,  in  particular,  the  U.S.  Army  Materiel  Command  (AMC) 
explosive  series,  “Properties  of  Explosives  of  Military  Interest”  [4],  Next,  we  describe  in  more 
detail  what  is  meant  by  the  term  “blast  wave.” 

When  energy  from  a  source  is  rapidly  released  into  the  surrounding  medium,  a  finite 
pressure  wave  of  compressible  gas  propagates  through  the  medium.  The  front  of  the  disturbance 
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steepens  as  it  propagates  forward.  The  air  is  said  to  be  shocked  up  as  it  passes  through  the  front, 
and  the  wave  is  termed  to  be  a  “shock  wave.”  The  shock  wave  front  propagates  at  supersonic 
speed  as  it  travels  through  the  undisturbed  air.  The  shock  wave  front  constitutes  a  jump 
discontinuity  in  the  properties  of  the  air  behind  the  front.  The  term  “blast  wave”  will  be  taken  to 
mean  a  shock  wave  generated  by  the  detonation  of  an  HE  charge  in  the  ambient  air  or  by  one  of 
the  energy  sources  previously  mentioned.  The  interest  in  our  present  study  is  confined  to  blast 
waves  generated  by  the  detonation  of  HE  charges. 

2.  Description  of  a  Generic  Blast  Wave 

2.1  Blast  Wave  in  Free  Air.  Let  us  consider  the  characteristics  of  an  ideal  blast  wave 
formed  in  the  air  by  some  energy  source.  Assume  that  an  HE  explosion  occurred  in  air,  that  it  is 
spherically  symmetric,  and  that  the  atmosphere  is  homogeneous.  It  follows  that  the 
characteristics  of  the  air  behind  the  shock  front  are  only  functions  of  the  distance  R  from  the 
source  of  the  disturbance  and  the  elapsed  time  t.  If  an  ideal  pressure  transducer  that  offers  no 
resistance  to  the  flow  behind  the  shock  front  was  located  a  short  distance  from  the  source  of  the 
HE  explosion,  then  this  transducer  will  record  the  variations  of  the  pressure  behind  the  shock 
wave  as  it  sweeps  over  the  transducer.  Tbe  record  produced  by  such  a  transducer  will  resemble 
that  shown  in  Figure  1.  For  some  time  after  the  explosion,  the  transducer  records  ambient 
pressure  po  until  the  shock  wave  arrival.  At  arrival  time  t„  the  pressure  rises  quite  abruptly 

(discontinuously  in  an  ideal  wave)  to  a  peak  value  p^  +  p^ .  The  pressure  then  decays  to 
ambient  in  a  total  time  t,  + 1”,  drops  to  a  partial  vacuum  of  amplitude  p^" ,  and  eventually  returns 

to  Po  in  total  time  tj,  +  T”  +  T.  The  quantity  p/  is  usually  termed  the  peak  side-on  overpressure 

or  merely  the  peak  overpressure.  Tlie  portion  of  the  time  history  above  the  initial  ambient 
pressure  po  is  called  the  positive  phase  of  duration  T*.  The  portion  below  po,  of  amplitude  p^ 

and  duration  T',  is  called  the  negative  phase.  Tlie  positive  and  negative  impulses  defined  by  die 
equations 
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Figure  1.  Ideal  Blast  Wave. 

C  =  wo  -Pol  <it  (1) 

and 

V  [Po-P(01dt  (2) 

are  important  blast  wave  parameters.  In  the  following,  we  shall  focus  our  attention  only  on  the 
characteristics  of  the  positive  phase  part  of  the  blast  wave  and  its  significant  parameters. 
Henceforth,  any  blast  wave  parameters  mentioned  will  refer  only  to  the  positive  phase  part  only. 
Li  order  to  describe  the  characteristics  of  an  ideal  blast  wave,  we  &st  need  a  functional 
description  of  its  pressure-time  history.  Several  functional  forms  have  been  proposed  for  that 
purpose.  A  two-parameter  form  approximates  the  wave  with  a  straight  line  that  preserves  two  of 
the  blast  wave  parameters  (e.g.,  p^  and  T*"),  as  is  the  case  in  equation  (3): 

p(t)  =Po  +p;(i-tAr*),  0<t<=T\  (3) 

where,  tj  is  assumed  to  be  zero  and  henceforth.  A  frequently  used  three-parameter  form  is  the 
“modified  Friedlander  equation”: 

p(t)  =  Po  +  p;  0  <  t  <=  T\  (4) 
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This  three-parameter  form  allows  one  to  match  any  three  of  the  four  blast  characteristic 
parameters,  P/,  T  %  ,  and  initial  decay  rate  dp/dt(t  =  0).  A  more  complex  four-parameter 

form  has  been  proposed  by  Ethridge  [5]  that  will  match  all  the  blast  parameters  previously 
shown.  Excellent  references  on  the  phenomenology  of  blast  waves  and  explosions  can  be  found 
in  AMCP  706-181  [6]  or  Baker  [7]. 

12  Shock  Wave  Reflection.  When  a  shock  wave  strikes  a  rigid  wall  or  a  solid  obstacle,  the 
properties  of  the  wave  are  changed  considerably  as  the  wave  reflects  from  the  obstacle  or 
diffracts  around  it.  Many  books  exist  that  discuss  the  subjects  of  wave  reflection  or  diffraction  in 
For  our  purposes,  we  shall  briefly  describe  the  phenomenon  of  normal  reflection  of  a 
plane  shock  wave  and  make  a  few  general  comments  about  shock  wave  reflections.  Figure  2 
illustrates  schematically  the  reflection  of  a  plane  incident  wave  I  from  a  rigid  wall.  The  shock 
front  moves  at  a  velocity  U  into  the  ambient  atmosphere  with  properties  po,  po,  0o,  and  Uq  =  0 
denoting  the  pressure,  density,  temperature,  and  particle  velocity.  The  flow  behind  the  shock 
front  is  characterized  by  the  pressure,  density,  temperature,  and  particle  velocity  denoted  by 
Po  -h  P„  Pj,  0„  and  u,.  At  the  rigid  wall,  the  particle  velocity  is  reduced  to  zero,  u,  =  0.  The 
reflected  wave  front  R  moves  away  from  the  wall  at  a  velocity  into  the  flow  field  behind  the 
incident  wave.  The  ovapressure,  density,  and  temperature  behind  the  reflected  wave  are 
substantially  increased  to  new  values,  Po  +  P,,  Pp  and  0,.  For  a  weak  shock  wave  (i.e.,  Pj «  po), 
the  reflected  overpressure  can  be  shown  to  be  P^  =  2  Pj,  and  the  wave  is  called  an  acoustic  wave. 
In  general,  the  reflected  overpressure  can  be  many  times  the  incident  overpressure.  For  an  ideal 
gas,  it  can  be  shown  that  P^  =  S-P^.  However,  it  was  shown  by  Doering  and  Burkhardt  [8]  and 
Shear  and  McCane  [9]  that  air  does  not  behave  as  an  ideal  gas  under  high-temperature  and 
high-pressure  conditions.  It  is  believed  that  under  such  conditions,  P,  can  be  as  large  as  20-P5. 
Shock  waves  can  also  reflect  obliquely  from  a  solid  obstacle  or  rigid  wall.  This  reflection 
phenomenon,  as  well  as  the  Mach  reflection  of  a  shock  wave,  are  highly  complex  and  are  beyond 
the  scope  of  this  study.  The  interested  reader  should  see  Baker  [7]  and  Courant  and  Friedrichs 
[10]. 
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Figure  2.  Normal  Reflection  of  a  Plane  Shock  Wave  From  a  Rigid  Wall. 

23  Blast  Parameters  Scaling  Laws.  Determination  of  the  blast  parameters  that 
characterize  the  blast  wave  from  an  explosive  charge  W  detonated  in  the  air  can  be  made  by 
experimental  means.  These  experiments  are  costly  and  time  consuming.  Computational 
methods  for  calculating  blast  wave  parameters  are  also  highly  complex,  time  consuming,  and 
expensive.  Early  on,  researchers  started  looking  for  scaling  methods  in  order  to  denve  from  then 
experiments  the  values  of  blast  parameters  for  many  more  similar  explosive  charges  of  different 
size  without  resorting  to  experiments  for  each  possible  variation  of  characteristic  parameters. 
Hopkinson  or  “cube  root”  scaling  is  most  widely  used  by  people  who  deal  with  blast  waves  from 
explosive  charges.  This  scaling  law  states  that  self-similar  blast  (shock)  waves  are  produced  in 
air  when  two  explosive  charges  of  the  same  geometry  and  composition  but  of  different  sizes  are 
detonated.  Symbolically  stated  as: 

p  =  p(Z),  T  =  T(Z),  U  =  U(Z),  C  =  C(Z) 

where, 

Z  =  R/E‘'3,  -c  =  and  C  =  • 


(5) 

(6) 
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p,  I,  U,  and  E  denote  the  pressure,  impulse,  velocity,  and  total  energy.  R  is  the  radial  distance 
measured  from  the  center  of  the  charge,  and  t  denotes  the  time.  Z,  t,  and  C  are  designated  scaled 
distance,  scaled  time,  and  scaled  impulse.  Hopkinson  scaling  has  been  verified  experimentally 
over  a  wide  range  of  scaled  distances  and  can  also  be  derived  from  dimensional  analysis 
considerations  [6].  Traditionally,  experimentalists  have  used  (instead  of  the  dimensionless 
scaled  variables  of  Z,  t,  and  C  defined  previously)  dimensional  scaled  variables  that  result  from 
the  previous  definition  by  replacing  the  variable  E  by  W,  the  charge  weight  in  those  definitions. 
Since  W  is  proportional  to  E,  Hopkinson  scaling  still  holds  in  this  form  that  is  most  commonly 
used.  This  scaling  states  that  pressures  and  velocities  are  identical  at  the  same  scaled  distances 
Z,  and  similarly  scaled  impulse  and  time  are  functions  of  Z  only.  This  scaling  law  is  commonly 
used  to  predict  the  blast  parameters  for  spherical  HE  charges  from  the  documented  properties  of 
Pentolite  by  Goodman  [2]  or  from  those  of  TNT  as  documented  by  Swisdak  [1].  Peak  shock 
overpressures,  impulse,  and  time  of  positive  durations  can  be  predicted  for  the  side-on  free  air 
values  or  for  the  normally  reflected  wave  parameter  values  from  established  standard  curves  and 
tables.  Hopkinson  scaling  is  also  used  to  predict  blast  parameter  values  for  cylindrical  HE 
charges  from  established  data  for  such  charges.  Predictions  of  blast  parameter  values  such  as 
overpressure  and  impulse  are  essential  for  predicting  the  blast  damage  to  the  various  structural 
components  of  a  target  subjected  to  HE  blast. 

In  order  to  account  for  the  variations  in  the  values  of  the  blast  parameters  with  altitude 
conditions,  Sachs  [1 1]  proposed  a  more  general  scaling  law  than  Hopkinson’s  scaling.  Sachs’ 
scaling  states  that  dimensionless  parameters  involving  pressure,  time,  and  impulse  can  be  defined 
in  combination  with  other  parameters  of  the  ambient  atmosphere  so  that  they  are  unique 
functions  of  a  dimensionless  distance  parameter.  Symbolically,  it  can  be  stated  that  the 
dimensionless  parameters 


p/p„,  I  .  V(E  .  Po')‘'’,  t  .  a,  .  p<;'’/E''’ 

are  unique  functions  of  the  dimensionless  scaled  distance; 


(7) 
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(8) 


1/3 /C  1/3 


R  .  Po‘'7E 


where  p,  I,  R,  and  t  represent  the  pressure,  impulse,  radial  distance  from  the  charge  center,  and 
the  time.  E  is  the  total  energy,  with  ao  and  po  being  the  sound  speed  and  the  ambient  atmospheric 
pressure.  Sperrazza  [12]  gave  a  derivation  of  Sachs  scaling  using  dimensional  analysis 
argument;  Dewey  and  Sperrazza  [13]  verified  the  scaling  experimentally.  Using  superscripts  on 
variables  to  denote  altitude  values  and  zero  subscripts  to  denote  ambient  atmosphere  values  of 
parameters,  Sachs  scaling  can  be  stated  by  the  following  equations: 

P^Vp(f^  =  p(°Vp(J°\  (9) 

ao^^  .  (10) 

and 

t(")  .  ao^^  .  {po^VE<^>}‘/3  ^  t(0)  _  ^(0)  _  {pW/E(0)}i/3.  (H) 

Equations  (9)-(l  1)  summarize  Sachs’  scaling  law.  These  equations,  together  with  equation  (8), 
can  be  used  to  predict  values  of  the  blast  parameters  at  altitude  condition  (h)  from  values 
measured  at  sea  level.  Values  of  blast  parameters  at  altitude  conditions  are  essential  for 
vulnerability/lethality  studies  of  aircraft  and  missiles  in  flight. 


3.  Illustrative  Experimental  Test  Results 


•  HE  Blast  vs.  Rolled  Homogeneous  Armor  (RHA)  Plates,  Aluminum  Plates,  Cylinders, 
Trucks,  Armored  Personnel  Carriers  (APCs),  Scaled  Chemical  Submunition  Cylinders, 
and  Aircraft.  Some  experimental  test  data  exist  on  the  damage  caused  by  HE  blast  waves  to 
various  structures.  Data  of  HE  blast  damage  to  RELA  plates  of  3/1 6-3/ 4-in  thickness  have  been 
documented  in  Kiwan  and  Goodman  [14].  Figures  3, 4,  and  5  show  the  various  types  of  damage 
observed  in  these  experiments.  The  data  gathered  in  these  tests  consisted  of  measurements  of 
deformation,  reduction  in  thickness,  failure  threshold  parameters  values,  and  hole  size 
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Figure  3.  Damage  to  a  3/16-in-Thick  RHA  Plate  by  a  5-Ib  Spherical  Pentolite  Charge  at  a  6.25-in  Standoff. 
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Figure  4.  Damage  to  a  3/8-in-Thick  RHA  Plate  by  a  2-lb  Spherical  Pentolite  Charge  at  a  2.5-m  Standoff. 


Figure  5.  Damage  to  a  3/8-m-Thick  RHA  Plate  by  a  5-Ib  Spherical  Pentolite  Charge  at  a  3.75-m  Standoff. 


measurements.  The  experimental  test  data  were  used  to  construct  a  mathematical  model  for  the 
response  of  armor  to  HE  blast  loading.  The  developed  model  will  be  discussed  in  section  4.2. 
Figure  3  shows  the  HE  blast  damage  of  an  RHA  plate  3/16  in  thick  from  the  detonation  of  a  5-lb 
Pentolite  charge  at  a  standoff  of  6.25  in.  The  observed  damage  consisted  of  deformation  and 
buckling.  Figure  4  shows  the  more  severe  damage  of  material  failure  with  three  cracks 
propagating  from  the  plate’s  center.  Figure  5  shows  a  more  severe  type  of  damage  of  a  hole 
opening  up,  with  cracks  propagating  toward  the  plate’s  boundary  and  petaling. 

Some  experimental  data  exist  on  the  response  of  aluminum  2024  plates  to  blast  loading  from 
2-lb  Pentolite  charges.  The  available  data  are  limited  to  plate  thicknesses  of  0.090  and  0.125  in. 
Figure  6  shows  the  blast  damage  from  a  2-lb  Pentolite  charge  detonated  at  a  standoff  of 
11  9/16  in  from  an  aluminum  2024  plate  0.090  in  thick.  The  plate  failed  in  shear  along  the  entire 
clamped,  supported  boundary.  Figure  7  shows  a  similar  experiment  but  at  a  standoff  of  22  in. 
Buckling  and  tearing  were  observed  at  the  central  bolts  in  addition  to  the  deformation  and  the 
tearing  at  the  comers.  Figure  8  shows  a  similar  test  for  a  plate  of  0.125-in  thickness  at  a  charge 
standoff  of  13  in.  The  damage  observed  in  these  experiments  differed  in  nature  from  the  damage 
observed  on  the  RHA  plates.  The  aluminum  2024  deformed  plastically  but  then  failed  in  shear  at 
the  boundary  starting  at  the  plate’s  comers.  The  RHA  plates  failed  in  tension  from  the  plate’s 
center.  Schuman  [15, 16]  conducted  extensive  experiments  on  the  blast  damage  to  aluminum 
and  steel  cylinders.  The  results  of  these  experiments  and  research  are  documented  in  noted 
references.  The  resulting  damage  model  of  Schuman  will  be  described  briefly  in  section  4.3. 
Baker  et  al.  [17]  investigated  the  response  of  cantilever  beams  to  air  blast  loading. 

Experimental  test  data  of  the  response  of  stmctures  of  various  weapon  systems  to  HE  blast 
loading  exist,  although  sometimes  not  in  documented  form.  HE  blast  tests  were  made  against 
U.S.  Army  tmcks,  APCs,  tank  tracks  and  bottoms,  fixed-wing  and  rotary-type  aircraft,  and 
various  missile  stmctural  components.  Figure  9  shows  preparation  for  the  test  of  blast  effects 
from  an  18.25-lb  cylindrical  Pentolite  charge  at  a  standoff  of  10  ft  from  the  skin  of  a  B-57 
aircraft.  Figure  10  shows  the  sustained  blast  damage,  and  Figure  1 1  shows  a  closeup  view  of  the 
test  damage.  Five  HE  blast-damage  tests  were  made  on  B-57  aircraft.  The  test  results  and  data 
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Figure  6a.  Pretest  Setup  of  a  0.090-in-Thick  Aluminum  2024  Plate  at  an  11 9/16-in 
Standoff  From  a  2-lb  Spherical  Pentolite  Charge. 


Figure  6b.  After-Test  Damage  to  a  0.090-in-Thick  Aluminum  2024  Plate  Due  to  the  Blast 
From  a  2-lb  Spherical  Pentolite  Charge  at  an  11  9/16-in  Standoff. 
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Figure  7.  Damage  to  a  0.125-in-Thick  Aluminum  2024  Plate  Due  to  the  Blast  From  a 
2-lb  Spherical  Pentolite  Charge  at  a  13-in  Standoff. 


Figure  8.  Damage  to  a  0.090-in-Thick  Aluminum  2024  Plate  Due  to  the  Blast  From  a  2-lb 
Spherical  Pentolite  Charge  at  a  22-in  Standoff. 


13 


14 


Figure  9.  Preparations  for  Blast  Test  of  an  18.25-lb  Bare  Pentolite  Cylinder  Charge  vs.  B-57  Aircraft  at  a  10-ft  Standoff. 


e  of  a  B-5' 
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analysis  are  shown  in  Mayerhofer,  Kiwan,  and  Williams  [18].  Figure  12  shows  the  damage 
sustained  by  an  OH-58  helicopter  due  to  the  HE  blast  from  an  8-lb  spherical  Pentolite  charge 
detonated  at  a  standoff  of  7  ft  4  in  from  the  skin  of  the  cargo  cabin  door.  Figure  13  shows  the 
damage  sustained  by  Hound  Dog  (AGM-28)  missile  engine  inlet  when  subjected  to  the  blast 
from  a  surface  burst  of  an  8-lb  spherical  Pentolite  charge  at  a  distance  of  7  ft  from  the  skin  of  the 
inlet.  The  data  from  similar  tests  are  in  Kiwan  [19].  Because  HE  blasts  are  more  effective 
against  light  structures,  extensive  experimental  blast  data  are  available  against  various  aircraft 
and  missile  structures.  Test  data  of  the  Redeye  missile  (XM221)  against  various  aircraft  and 
missile  structures  can  be  found  in  the  U.S.  Army  Test  and  Evaluation  Command  firing  record 
R-3708  [20].  Damage  data  to  the  UH-1  helicopter  from  the  blast  of  small  C-3  charges  are  also 
available,  as  well  as  the  blast  damage  data  on  a  number  of  other  aircraft.  Finally,  Figure  14 
shows  the  damage  sustained  by  scaled  steel  cylinders  simulating  chemical  weapons  containers 
due  to  the  blast  from  a  20-lb  cylindrical  Pentolite  charge. 

4.  HE  Blast  Damage  Models 

4.1  Blast  Hole  Model.  Wolk  and  Wilner  [21]  developed  a  formula  for  the  blast  hole  radius 
DR,  in  an  armor  plate  of  thickness  T,  due  to  the  blast  from  an  HE  charge  W,  detonated  at  a 
standoff  distance  S  from  the  plate.  They  assumed  that  an  RHA  plate  of  thickness  T  will  fail 
when  subjected  to  a  critical  overpressure  P^.  Furthermore,  they  assumed  that 

P,  =  a.T\  (12) 

where  a  and  b  are  constant  parameters.  By  fitting  the  side-on  overpressure  data  as  a  function  of 
scaled  distance,  they  derived  an  equation  of  the  form 

p=k.Z‘‘,  (13) 
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Figure  12a.  Preparation  for  an  HE  Blast  Test  of  an  8-Ib  Spherical  Pentolite  Charge  vs.  an 
OH-58  Helicopter  at  a  7-ft,  4-in  Standoff. 
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Figure  13.  Damage  to  a  Hound  Dog  Missile  (AGM-28)  Engine  Inlet  Due  to  the  Blast  From  an  8-lb  Surface  Burst  at  a  7-ft 
Standoff. 


where  k  and  d  are  constant  parameters  and 


Z  =  (14) 

The  blast  from  the  charge  W  will  generate  a  hole  in  the  plate  of  radius  DR,  where  the  normal 
component  of  the  blast  overpressure  is  equal  to  the  critical  overpressure  for  that  plate 
thickness  T.  By  combining  equations  (12)  and  (13),  they  arrived  at  an  equation  for  the  blast  hole 
radius  DR  as  a  function  of  W,  S,  and  T.  The  formula  can  be  expressed  in  functional  form  as: 


DR  =  F(W,  S,T). 


(15) 


Comparison  of  predictions  from  equation  (15)  with  test  results  of  holes  generated  by  the  blast 
from  large  bombs  showed  qualitative  agreements.  This  model  has  been  used  in  some  ship 
vulnerability  prediction  codes. 


4.2  RHA  Plates  Blast  Damage  Model.  Kiwan  used  the  experimental  RHA  test  data  for  HE 
blast  damage  discussed  in  the  previous  section  and  documented  in  Kiwan  and  Goodman  [14]  in 
combination  with  similarity  and  dimensional  analysis  to  derive  a  mathematical  model  for  HE 
blast  damage  to  RHA  plates  [22].  The  analysis  was  later  extended  to  apply  to  moving  HE 
charges  and  documented  in  Kiwan  [23].  The  derived  mathematical  equations  are 

d  =  a .  .  S^“*),  (16) 


S  ■^5521/-p0il49 


(17) 


and 


DR  =  [A^ .  .  'r‘-«233)  -S2]°•^ 


(18) 


where  d  denotes  the  maximum  plate  deformation,  W,  S,  and  T  are  the  charge  weight,  the  standoff 
of  the  charge  center  from  the  plate,  and  the  plate  thickness,  a  and  A  are  constant  parameters 
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whose  values  depend  on  the  system  of  units  used.  denotes  the  critical  standoff  for  charge  W, 
at  which  the  plate  fails  and  cracks  appear. 

a  =  27392697  and  A  =  1.4478  x  10“,  (19) 

when  the  millimeter  is  used  for  the  unit  of  length  and  the  kilogram  is  used  for  the  unit  of  mass. 

a  =  0.8631  and  A  =  2.0317,  (20) 

when  the  inch  and  pound  are  used  as  units. 

Predictions  made  with  this  model,  as  expressed  with  equations  (16)-(18),  compare  well  with 
the  experimental  data  within  the  range  of  its  validity,  i.e., 

0.1875  <=T<=  0.75  in,  (21) 

and 

2.0  <-  W  <=  8.0  lb.  (22) 

For  more  details  on  this  model  and  comparison  of  predictions  with  experimental  results,  we  refer 
the  interested  reader  to  Kiwan  [22, 23]. 

4,3  Cylindrical  Shells  Blast  Damage  Model.  Schuman  conducted  extensive  experiments 
on  the  response  of  aluminum  and  steel  cylindrical  shells  to  blast  loading.  TTie  data  from  these 
experiments  were  used  to  formulate  a  model  for  the  damage  to  such  shells  [15, 16].  For  lateral 
loading  of  cylinders  by  the  blast  wave,  a  critical  overpressure  was  determined  for  a  shell  to 

sustain  a  5%-10%  permanent  deformation  of  its  diameter.  The  required  critical  overpressure 
was  found  to  be  a  function  of  the  shell  material,  length-to-diameter  ratio,  shell  diameter, 
thickness,  and  the  charge  weight.  This  functional  dependence  can  be  expressed  by  the  equation: 

Pcr  =  FuD-FD-FT-Fw»  (23) 
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where  Fl/j,  .  .  F^.  .  and  F^  are  the  factors  for  length-to-diameter  ratio,  the  diameter,  the  shell 

thickness,  and  the  charge  weight.  For  a  given  case,  these  factors  are  determined  from  graphs  that 
are  given  in  Schuman  [15].  For  longitudinal  loading  of  these  shells,  the  required  critical 
overpressure  P„  can  be  expressed  as: 


Pcr  = 


g-P 


cr’ 


(24) 


where  g  is  a  multiplication  factor: 


(25) 

(26) 

This  model  has  been  used  in  missile  vulnerability  studies. 


and 


g  =  6.0  for  steel 


g  =  2.0  for  aluminum. 


5.  HE  Blast  Damage  Scaling 


5.1  The  Johnson  Relationship.  Johnson  [24]  developed  a  relation  between  the  explosive 
HE  charge  weight  W,  and  the  range  R,^,  at  which  a  certain  level  of  target  response  or  damage  due 
to  the  HE  blast  wave  is  achieved.  Johnson  examined  the  blast  damage  data  for  various  types  of 
targets  (structures)  and  made  the  following  observations.  For  each  combination  of  charge  weight 
W  and  distance  R,  one  can  associate  values  of  pressure  P  and  impulse  I,  and  conversely.  If  one 
plots,  for  a  given  target  or  structure,  the  set  of  P  and  I  data  that  will  result  in  an  equivalent  level 
of  damage,  then  such  data  fall  on  a  curve  similar  to  the  one  in  Figure  15,  which  can  be  assumed 
to  approximate  a  rectangular  hyperbola.  The  target  set  that  Johnson  examined  consisted  of  wire 
drag  gauges,  aluminum  beams,  aluminum  cylinders,  dish  radar  antenna,  B-29  fuselage  section, 
and  2-1/2-ton  trucks.  Johnson  suspected  that  for  a  given  target  and  level  of  response  a 
relationship  should  exist  between  the  explosive  HE  weight  W,  and  the  range  R^  at  which  that 
level  of  response  is  achieved.  Hence,  by  choosing  an  arbitrary  but  constant  explosive  weight  and 
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Figure  15.  An  Isodamage  (Pressure-Impulse)  Curve. 


referencing  all  data  to  it,  we  should  arrive  at  a  constant  range  ratio  for  that  weight  W.  For  the 
constant  weight,  Johnson  chose  W  =  100  lb  and  defined 


For  the  various  considered  targets,  Johnson  made  a  least-squares  fit  of  C,,  vs.  W,  to  the 
functional  form 


-  a .  W^  (28) 

Figure  16  shows  the  least-squares  fit  of  the  data  used  by  Johnson.  The  a  and  b  were  determined 
to  be  7.64  and  -0.435.  Equation  (28)  becomes; 


C^  =  7.64W“"‘^^ 


(29) 
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explosive  weight,  w  lbs. 

Figure  16.  A  Least-Squares  Fit  of  vs.  Explosive  Weight  for  an  Arbitrary  Response 
Level. 

Let  Rj,  Rj  be  the  standoif  distances  corresponding  to  charge  weights  Wj  and  Wj.  Equation  (29) 
then  gives  the  scaling  relationship: 

Rl/R2  =  (W,/W/•^3^  (30) 

Equation  (30)  can  be  used  to  scale  from  one  vulnerability  data  point  corresponding  to  charge 
weight  Wj  to  other  charge  weights.  Equation  (30)  has  been  generalized  to: 

R,/R2  =  (WiAV2)'.  (31) 

The  value  of  the  parameter  c  in  equation  (3 1)  depends  on  the  range  of  values  of  W,  and  Wj. 
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c  -  0.333 

for 

O 

11 

o 

for 

c  -  0.435 

for 

c  =  0.49 

for 

c-0.6 

for 

Wj  in  the  kiloton  range 
500  lb  <  Wj  <  kiloton 
101b<Wi<5001b 
1  lb  <  Wj  <  10  lb 
W;  <  1  lb 


(32) 


5.2  Isodamage  Curve  Method.  The  concept  of  an  isodamage  curve  expressed  in  terms  of 
the  parameters  of  peak  oveq)ressure  P  and  impulse  I  was  first  proposed  by  Speirazza  [25].  This 
concept  is  based  on  the  experimental  observations  that  some  parameters  of  the  response  of 
simple  stmctures,  such  as  cantilever  beams,  to  transient  loads  such  as  maximum  displacement, 
are  a  function  of  the  peak  load  and  impulse  delivered  to  the  target.  Plots  of  the  overpressure  P 
vs.  the  impulse  I  data  for  equal  values  of  maximum  deflection  fall  on  an  isodamage  curve  similar 
to  the  one  shown  in  Figure  15.  Therefore,  one  can  assume  that  such  a  relationship  imght  hold  for 
more  severe  damage  of  complex  structures.  It  is  possible  to  estimate  for  a  given  explosive 
weight  W,  a  combination  of  overpressure  P,  and  impulse  I  that  will  result  in  a  specific  level  of 
damage  to  a  target.  The  P  -  I  curve  shown  in  Figure  15  separates  the  plane  into  two  distinct 
regions.  Data  points  to  the  right  of  the  curve  represent  loading  that  will  cause  more  severe 
damage.  Points  to  the  left  of  the  curve  represent  loading  that  will  cause  less  severe  damage  to  the 
structure.  Greenspon  [26]  gave  a  theoretical  analysis  basis  for  the  existence  of  isodamage  P  -  I 
curves.  The  curve  can  be  assumed  to  be  hyperbolic  and  can  be  represented  by  the  equation. 


(P-Pcrit.)a-Icrit.)  =  C. 


(33) 


C  is  a  constant  parameter  whose  value  depends  on  the  target’s  characteristics.  Pen,,  and  lerf,,  are 
two  parameters  representing  the  curve  asymptotes.  Determination  of  a  specific  isodamage  curve 
requires  three  data  points.  However  using  the  scaling  relationship  (equation  31)  enables  one  to 
determine  such  a  curve  once  we  have  a  single  data  point.  The  P  -  I  curve  can  be  used  to  scale  to 
other  overpressure  P  and  impulse  I  values  that  will  result  in  the  same  level  of  damage  to  a  given 

target. 
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5.3  Other  Scaling  Relationships.  Equation  (17)  gives  the  critical  standoff  S^,  for  an  HE 
charge  of  weight  W  from  an  RHA  plate  of  thickness  T  in  order  for  the  plate  to  fail  and  crack  due 
to  the  HE  blast  loading.  One  might  suspect  that  this  relation  applies  to  other  stractures.  Hence, 
dropping  the  subscript  c  and  using  numerical  subscripts  to  describe  parameters  associated  with 
various  data  sets,  one  deduces  from  equation  (17)  the  following  relations: 

For  T  =  constant,  S/Sj  =  .  (34) 

Equation  (34)  is  a  scaling  relation  analogous  to  equations  (30)  or  (31).  Equation  (30)  is  based  on 
a  least-squares  fit  of  the  data  set  used  by  Johnson  and  documented  in  his  BRL  report  [24],  while 
equation  (34)  is  based  on  the  RHA  plates  data  set  documented  in  Kiwan  and  Goodman  [14]. 
Equation  (34)  can  be  used  to  scale  for  a  given  target  from  one  charge  weight  to  another. 

Similarly,  we  can  deduce  from  equation  (17)  the  following  relation: 

For  W  =  constant,  S/S^  =  (T/rj)®-^!^^  (35) 

Equation  (35)  can  be  used  to  scale  from  one  target  to  another  for  the  same  constant  charge 
weight  W.  Equations  (34)  and  (35)  have  been  used  to  scale  experimental  results  of  blast  damage 
of  steel  cylinders  to  full-size  submunition  containers. 

6.  HE  Blast  Vulnerability  Models 

6.1  The  Blast  Correction  Factor  Method.  R.  Kirby  and  H.  Ege  attempted  to  account  for 
the  contribution  made  by  HE  blast  loading,  case  fragments  perforations,  and  transmitted  shock 
from  impacting  explosive  munitions  to  the  vulnerability  of  lightly  armored  vehicles.  Based  on 
some  combat  data,  they  estimated  the  probability  of  kill  due  to  these  damage  mechanisms  as  a 
function  of  the  explosive  charge  weight  of  the  munition  and  the  mean  armor  thickness  of  the 
armored  vehicle.  This  probability  of  kill  denoted  by  Pbcs  is  shown  in  Figure  17.  This  probability 


27 


Figure  17.  Probability  of  Kill,  Due  to  Blast,  Case  Fragments,  and  Shock. 
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of  kill  is  combined  with  the  probabilities  of  kill  for  other  damage  mechanisms  present,  such  as 
the  probability  of  kill  due  to  the  shaped  charge  jet  through  the  survival  rule.  Thus  one  arrives 
at  a  total  probability  of  kill  P,o,  for  the  munition. 


P,,,  =  1  -  (1  -  PJ  (1  -  Pi,,) 


(36) 


This  method  of  accounting  for  the  contribution  of  blast  effects  to  the  vulnerability  of  armored 
vehicles  provides  quick  vulnerability  estimates  to  threat  munitions.  The  main  disadvantage  of 
this  method  is  the  small  amount  of  test  data  on  which  it  is  based. 

6.2  The  Critical  Impulse  in  Critical  Time  Method.  Robert  Sewel  of  the  U.S.  Naval 
Ordnance  Test  Station  (NOTS)  formulated  “A  Blast  Damage  Criterion”  and  a  vulnerability 
model  for  structures  due  to  HE  blast  effects  [27].  According  to  Sewel’ s  model,  a  structure 
sustains  permanent  damage  when  a  critical  impulse  is  delivered  to  the  structure  in  a  critical  time 
period.  A  formula  is  given  for  computing  the  critical  impulse  I,.; 

I.  =  .  p  .  T,  (37) 

where  1,,  V^,  p,  and  T  are  the  critical  impulse,  critical  impact  velocity,  and  material  density  of  the 
target  and  thickness,  respectively.  The  critical  impact  velocity  is  related  to  the  critical  particle 
velocity  in  the  material,  which  is  that  relative  velocity  between  two  particles  in  the  material,  and, 
if  exceeded,  causes  the  particles  to  behave  as  independent  particles  and  no  longer  as  part  of  the 
whole  system.  Furthermore,  it  was  assumed  that  the  critical  time  period  during  which  the  critical 
impulse  needs  to  be  delivered  to  the  target  to  achieve  the  desired  level  of  damage  is  pe/4,  where 
pe  is  the  period  of  the  fundamental  resonant  frequency  of  the  structure. 

The  application  of  this  damage  criterion  in  vulnerability  studies  was  achieved  by  computing 
the  side-on  and  the  normally  reflected  impulse  for  the  HE  blast  wave  and  assuming  that: 

15"“  P.  dt  >  ■  for  K,  -  1  (38) 
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and 


/r  P,  dt  <  =  I,  for  Kp  =  0, 


(39) 


where  P,  =  side-on  pressure,  psi,  Pj «  reflected  pressure,  and  psi  and  pe/4  are  the  quarter  period  in 
milliseconds.  Whenever  pe/4  >  'F,  the  integration  is  taken  to  T*,  the  time  of  positive  duration  of 
the  blast  wave,  t  =  time.  Kp  =  kill  probability  due  to  HE  blast.  Equations  (38)  and  (39)  define 
the  100%  and  0%  kiU-level  probabilities.  Once  the  radii  where  the  side-on  and  reflected  impulse 
equal  the  critical  impulse  1^  are  determined,  these  radii  are  defined  as  the  100%  and  0%  kill 
probabilities.  A  straight  line  is  drawn  between  those  two  data  points  in  the  (R,  Kp)  plane,  where 
R  is  the  radial  distance.  Sewel  used  this  model  and  computed  the  probabilities  of  kill  for  a 
number  of  military  systems  due  to  the  blast  from  large  HE  bombs.  Computations  were  made  for 
a  parked  aircraft,  a  gun  laying  radar  target,  a  missile  control  radar,  a  tmck,  an  APC,  and  a  tank. 
The  bombs  used  in  the  calculation  were  the  250-lb  MK-81,  the  500-lb  MK-82,  and  the  1000-lb 
MK-83,  and  the  U.S.  Air  Force  (USAF)  750-lb  T-54  bomb.  Experimental  values  of  parameters 
were  used  in  the  calculations  when  available.  Values  for  the  remaining  parameters  were  obtained 
from  tables  or  theory.  The  critical  impact  velocity  was  assumed  to  be  200  ft/s  for  steel  and 
300  ft/s  for  aluminum.  Here  we  present  some  of  the  results  of  these  calculations  as  reported  by 
Sewell  [27].  Figures  18-21  show  the  results  of  the  computation  for  four  of  these  targets. 

Application  of  the  Sewel  model  proved  to  be  complicated.  Computation  of  the  critical 
impulse  and  time  for  a  given  target  proved  to  be  a  difficult  task.  Similarly,  the  estimation  of  the 
critical  particle  velocity  for  a  given  material  is  not  an  easy  task.  These  parameters  often  had  to  be 
determined  from  experiments  that  are  rather  costly  and  defeat  the  purpose  of  the  method. 

6,3  Blast  Kill  Radius  Method.  Researchers  at  the  Naval  Surface  Warfare  Center  (NSWC) 
formulated  an  empirical  method  for  calculating  blast  kill  radii  for  missiles.  The  method  is  based 
on  test  data  to  determine  the  blast-critical  impulse  threshold  for  obtaining  catastrophic  structural 
kill  of  a  missile  and  applying  free-air  blast-measured  relationships  for  Pentolite  charges  to 
determine  where  the  required  critical  impulse  could  be  obtained.  A  brief  description  of  the 
method,  assumptions  made,  and  the  computational  steps  will  be  given  as  follows. 
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Figure  18.  Kill  Probability  vs.  Distance  for  Parked  Aircraft  Target 


Figure  19.  Kill  Probability  vs.  Distance  for  Gun-Laying  Radar  Target 
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Figure  20.  Kill  Probability  vs.  Distance  for  Truck  Target. 


Figure  21.  Kill  Probability  vs.  Distance  for  Tank  Target 
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Based  on  the  review  of  test  data  [28, 29,  30],  it  was  assumed  that  the  critical  damage  required 
to  provide  a  structural  kill  of  a  missile  and  cause  loss  of  control  and  ground  impact  within  30  s  is 
rupture  of  the  missile’s  skin.  Based  on  test  data  for  air  targets  with  skin  thicknesses  in  the  range 
of  0.040-0.075  in,  it  was  determined  that  a  reflected  impulse  of  100-180  psi-m/s  is  required  to 
cause  skin  rupture  in  this  range  of  thicknesses.  A  fit  of  the  experimental  test  data  of  critical 
impulse  vs.  skin  thickness  for  targets  with  aluminum  skin  gave  the  following  relationship: 

I,  =  2709.  T  psi-m/s.  (40) 

Equation  (40)  is  approximately  correct.  T  is  the  skin  thickness  in  inches.  In  order  to  use  this 
equation  in  vulnerability  assessments  to  a  threat,  one  has  to  determine  first  the  equivalent  bare 
charge  weight,  W^,  of  the  threat.  can  be  determined  from  the  modified  Fanno  formula; 

Wg  =  Explosive  Weight  [0.6  +  0.4  (1+2  Case  Weight/Explosive  Weight)"^].  (41) 

Next,  one  calculates  the  scaled  critical  impulse  of  the  target  and  equates  it  to  the  scaled  reflected 
impulse. 


2709  .  (42) 

From  the  curves  for  scaled  impulse  data  for  the  geometry  of  the  threat  charge,  one  determines  the 
scaled  distance  at  which  this  impulse  is  achieved.  The  lethal  blast  kill  radius  is  then  defined  as; 

R  =  (R/W^^)  (W,*'^).  (43) 

This  method  can  be  used  to  provide  quick  approximate  estimates  of  blast  kill  radii  for  aircraft 
and  missiles  in  the  absence  of  test  data. 

6.4  The  Lethal  Blast  Contours  Method.  Historically,  the  vulnerability  of  an  aircraft  to  the 
blast  from  an  HE  projectile  was  represented  by  defining  a  three-dimensional  volume  around  the 
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aircraft,  such  that  a  kill  of  the  aircraft  occurs  if  the  HE  charge  detonates  inside  that  volume. 
Visualization  of  this  volume  (termed  the  lethal  blast  volume)  was  achieved  by  taking  plane  cross 
sections  through  the  aircraft  and  the  surrounding  volume.  The  contour  created  by  the  intersection 
of  the  plane  and  the  lethal  volume  is  termed  the  lethal  blast  contour.  Mayerhofer  (formerly  of 
BRL)  formalized  this  process  for  aircraft  by  constructing  such  contours  in  four  orthogonal 
planes.  Three  of  the  planes  are  mutually  orthogonal.  TTiese  planes  are  usually  designated  as 
A,  B,  C,  and  D.  Plane  D  contains  the  aircraft  main  axis  and  is  parallel  to  the  wing  span.  Plane  C 
also  contains  the  aircraft  main  axis  and  is  perpendicular  to  plane  D.  Plane  A  is  perpendicular  to 
the  aircraft  main  axis,  but  intersects  some  critical  components  or  subsystems  in  the  forward 
section  of  the  aircraft  such  as  the  crew  compartment  or  the  wing  span.  Plane  B  is  parallel  to 
plane  A  but  intersects  some  critical  components  or  subsystems  in  the  rear  of  the  aircraft  such  as 
the  horizontal  or  vertical  stabilizer. 

The  computation  of  the  lethal  blast  contours  in  these  coordinate  planes  A,  B,  C,  and  D  for  an 
HE  charge  W  is  accomplished  as  follows:  For  the  kill  category  under  consideration,  an  analysis 
is  conducted  that  identifies  the  aircraft  critical  components  for  that  category  that  can  be  killed  by 
the  external  blast.  For  HE-blast  kills,  the  aircraft  structure  is  considered  a  critical  component. 
Next,  for  each  of  the  critical  components,  a  determination  is  made  of  the  type  and  extent  of 
damage  that  will  produce  a  complete  dysfunction  of  that  component.  For  the  plane  under 
consideration  and  each  of  the  critical  components  that  can  be  damaged  by  the  detonation  of  the 
HE  charge  W,  in  the  plane  or  a  small  neighborhood  of  it,  one  determines  the  maximum  radial 
distance  from  the  component  at  which  such  a  detonation  results  in  the  desired  level  of  damage. 
This  radial  distance  is  defined  as  the  lethal  blast  radius  R  for  that  component.  At  each  of  these 
critical  components,  a  circular  arc  of  radius  R  is  drawn.  The  lethal  blast  contour  in  the  plane  is 
obtained  by  drawing  the  envelop  of  these  arcs.  The  lethal  blast  volume  is  the  envelop  generated 
by  all  such  contours.  Figures  22-24  show  the  lethal  blast  contours  for  the  AH- IQ  helicopter  and 
for  a  series  of  TNT  charges  drawn  in  planes  A,  B,  C,  and  D. 
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Fieare  22  AH-IQ  Attritton  Kill  External  Blast  Contours  at  Sea^Uvel  tor  Equivalent 
Figure  22.  ^  3,.,^  ^  p,3„e. 


7.  Conclusions 

The  previously  described  methods  give  a  brief  overview  of  HE  blast  phenomenology,  HE 
blast  damage  mechanisms,  the  vulnerability  of  stmctutes  to  HE  blast,  and  vulnerabihty 
assessment  methods.  Hie  subjects  of  material  faUure  and  damage  assessments  are  highly 
complex  and  are  poorly  understood.  A  considerable  amount  of  research  of  material  failures  to 
various  damage  mechanisms  is  being  conducted  especially  on  composite  materials. 
Unformnately,  very  little  support  is  being  provided  for  studying  mamrial  faUure  and  damage  due 
,0  HE  blast  damage  mechanisms.  This  report  gives  a  brief  description  of  empirical  and 
engineering  methods  of  HE  blast  damage  and  assessment  methods  but  does  not  descnbe  any 
computer  codes  that  exist  on  the  subject  However,  this  report  states  that  the  accuracy  of 
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Figure  23.  AH-IQ  Attrition  Kill  External  Blast  Contours  at  Sea  Level  for  Equivalent  Charge  Weights  2.5-,  5-,  10-,  20-,  and 
30-lh  TNT;  C  Plane. 
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Figure  24.  AH-IQ 
30-lb  T] 


continuum  mechanics  codes  in  simulating  such  problems  relies  upon  the  material  failure  criterion 
contained  in  the  code  and  how  closely  it  simulates  material  behavior  under  these  loading 
conditions.  The  aforementioned  account  is  by  no  means  inclusive  of  all  available  models  and 
methods.  The  author  mentions  here,  in  passing,  a  model  for  material  failure  of  tank  bottoms  due 
to  a  mine  blast  formulated  by  Haskell  [31]  that  is  based  on  energy  considerations. 
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List  of  Symbols 


T,C 

d 

DR 

E 

^L/D’  Pd»  Pw 
Ic.V, 

C.  (i;) 

Kp 

Po>  Po»  ®0>  ^ 
Pbcs 
PcP. 
pe 

Pp  Pr»  ®r’ 

p;.  (p;) 

Ps»  Ps’  ®s’ 

Psc 

p(t) 

p 

^  tot 

R 


Scaled  time  and  scaled  impulse. 

Maximum  plate  deformation. 

Blast  hole  radius. 

Total  energy. 

Factors  of  length-to-diameter  ratio,  diameter,  thickness,  and  weight. 

Critical  impulse  and  critical  impact  velocity. 

Overpressure  and  ambient  atmosphere  sound  speed  at  altitude  h. 

Positive  (negative)  side-on  impulse. 

Kill  probability. 

Ambient  atmospheric  pressure,  density,  temperature,  and  sound  speed. 
Probability  of  kill  due  to  blast,  case  fragments,  and  shock. 

Critical  overpressure. 

Period  of  the  natural  resonant  frequency  for  the  structure. 

Overpressure,  density,  temperature,  and  particle  velocity  behind  the  reflected 
wave  front. 

Peak  side-on  positive  (negative)  phase  overpressure. 

Free  air  side-on  pressure,  density,  temperature,  and  particle  velocity. 
Probability  of  kill  due  to  shaped  charge  jet. 

Overpressure  at  time  t. 

Total  probability  of  kill. 

Radial  distance. 
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s,  s, 

T 

t, t, 

T^  (T-) 
W 

z 


Standoff  distance  and  critical  standoff  distance. 
Plate  or  skin  target  thickness. 

Time  variable  and  time  of  arrival  of  the  wave. 
Time  of  positive  (negative)  duration  of  the  wave. 
Explosive  weight. 

Scaled  distance. 
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